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I .  OBJECTIVE 

A  full  scale  fire  model  for  an  enclosed  environment 
needs  a  quantitative  description  of  the  fire  source  itself. 
This  is  because  of  the  time  dependent  power  output  of  a 
flame  arising  from  oxygen  depletion. 

The  objective  of  this  research  was  to  develop  a 
quantitative  understanding  of  the  behavior  of  simulated  gas 
flames  and  liquid  pool  fires  of  low  initial  momentum 
(buoyancy  controlled)  in  a  two-dimensional,  axi-symmetric 
situation.  Specifically,  we  predicted  velocity,  temperature 
and  species  concentration  fields  as  a  function  of  axial  and 
radial  coordinates.  Knowing  these  fields,  we  can  estimate 
the  quantity  of  feed  back  energy  which  in  turn  determines 
the  mass  burning  rate  of  a  liquid  pool.  This  enables  us  to 
determine  the  power  output  of  the  fire  source  in  question, 
as  a  function  of  the  chamber  environment. 
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II.  MODEL  FORMULATION 

A  brief  description  of  our  model  is  presented  here. 

The  model  to  be  described  has  been  developed  over  the  last  3 
years.  Those  features  incorporated  during  this  past  year 
are  indicated  in  Section  III  below. 

1 .Mathematical 

The  conservation  expressions  for  all  of  the  variables 
are  written  in  the  parabolic  form,  for  an  axi-symmetric 
flame  : 

Mass  : 

3<PU)  1  ^(ypV) 

-  +  -  -  =  0.  (I) 

^  z  y  ^  y 

Momentum  and  all  scalar  quantities  : 

b  $  1  ^  (y  (it  ^  ) 

pU  +  pV  --  =  -  --  -  --  +  R<?.  (2) 

dz  y^y(CT^t®y) 

The  variables  4> ,  the  Prandtl/Schmidt  numbers  cr<f,  and  the 
source  terms  have  been  tabulated  previously.  The 
effective  turbulent  viscosity  is  given  by 

dt  :  p  Cj  c!/t.  (3) 

The  boundary  conditions  irrespective  of  the  initial 
conditions  employed  are  : 

—  =  0  at  y=0  and  $  =  0  at  y -<&• 

by 

The  differential  equations  for  fuel  mass  fraction  (mfu) 
and  the  conserved  scalar  (mixture  fraction),  f,  are  solved 
within  this  framework.  The  mass  fractions  of  the  species  Os 
(1D02)  ,  CO2  (mCoi),  H20  (mH2o)  and  N2  (mn)  are  obtained  from 


the  auxiliary  relations: 


(<p  -  <po)/(<pt  -  <p 0)  ; 

(4. a) 

m  f  u  -  mo  2  / r  f  u 

(  4  .  b ) 

m f  u  +  mco2  /  ico2 

<4.c) 

m  f  u  +  mH2o  /  i b  2  0 

<4.d) 

mn2  ; 

<4.c) 

m  i  =  1 . 

(  4  .  d ) 

all  j 

Here,  the  subscripts  F  and  0  refer  to  the  fuel  bearing  and 
oxygen  bearing  streams  respectively;  and  the  i’s  are  the 
masses  of  the  species  CO2  and  H2O  produced  by  the  reaction 
of  a  unit  mass  of  fuel. 

The  specific  enthalpy,  h,  is  defined  by 

h  =  CP  T  +  (I-Xr)AH  m  f  u  (5) 

where  Cp  is  the  specific  heat,  AH  is  the  heat  of  combustion, 
and  Xr  is  the  radiative  fraction. 

The  temperature  dependence  of  the  laminar  viscosity,  u, 
is  given  by  the  relation, 

U  -  mo  (  T  /  Ta  )  0 • *  I  mi  ,  (6) 

where  Vo  is  the  viscosity  at  a  reference  temperature,  Ta. 

The  main  flow  field  is  assumed  to  be  parabolic.  This 
assumption  may  be  invalid  close  to  the  burner  exit,  hence 
some  error  may  exist  in  the  near  field  predictions. 

However,  similar  parabolic  flow  field  approximations  have 
been  made  by  other  authors.  The  usual  two  equation  model  is 
employed  to  describe  the  turbulence,  which  can  be 


4 


T" 
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characterized  by  the  turbulent  kinetic  energy,,  <c ,  and  the 
kinetic  energy  dissipation  rate,  e. 

2.  The  Combustion  Model 

The  chemical  kinetics  are  assumed  to  be  infinitely 
fast.  However,  the  microscale  processes  responsible  for  the 
turbulence  dissipation  and  molecular  mixing  in  the  turbulent 
field  are  typically  intermittent  and  certainly  not 
infinitely  fast. 

Magnussen  has  proposed  a  physical  and  conceptual  model 
for  the  interaction  of  the  turbulence  and  chemistry  for  both 
premixed  and  diffusion  flames.  In  this  model  the  chemical 
reactions  take  place  in  various  reactive  zones  or  fine 
structures.  The  reaction  rate  is  then  expressed  in  terms  of 
the  mass  transfer  into  the  fine  structures  from  the 
surrounding  fluid.  Thus  the  overall  reaction  rate  can  given 
by  the  expression, 

m  x 

Rf  a  :  - Cain  j  (7) 

(l-x  T*) 

where  Cain  is  the  smaller  of  Cfu  and  Co:/rf U(  which  indicate 
the  concentration  of  fuel  and  oxygen  in  the  flu'd 
respectively.  Here  rru  is  the  stoichiometric  oxygen 
requirement.  x  is  a  factor  between  zero  and  one  which  was 
introduced  by  Magnussen  for  premixed  flames  to  require  the 
presence  of  products  in  a  region  before  reaction  could 
occur.  Since  we  are  working  with  diffusion  flames  where 
this  is  inappropriate,  we  set  x  equal  to  one  . 
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The  quantity  m  is  the  mass  transfer  rate  between  the 
fine  structures  and  the  bulk  of  the  fluid;  they  can  be 
defined  in  terms  of  «c  and  c  and  the  kinematic  isotropic 
viscosity  v  as  follows; 


vc  1/4  e 

m  =  23.6  ( -  )  -  .  (8) 

K1  K 


T*  is  the  fraction  of  the  fluid  contained  within  the 
fine  structures  as  given  by 

ve  3  /  4 

T*  =  9.7  ( -  )  .  (9) 

K  * 

The  chemical  source  term  in  the  fuel  conservation  expression 
is  thus  given  by  Eq.(7).  This  chemistry-turbulence  source 
is  attractive  both  in  terms  of  its  applicability  to  fire 
problems,  and  to  both  premixed  and  diffusion  flame 
situations . 

3.  Radiation 

The  lack  of  understanding  and  development  of  radiative 
heat  transfer  models  in  fires  is  primarily  due  to  the  lack 
of  accurate  flame  property  data  needed  as  primary  input 
parameters . 

Considering  the  self  absorbing  and  emitting  non- 
homogeneous  media  in  fires,  as  well  as  in  the  present  pool¬ 
like  gas  flames,  we  chose  to  use  a  flux  model  for  radiation 
modeling.  Although  this  model  has  been  utilized  for  applied 
combustion  problems  such  as  in  cylindrical  furnaces,  the 
applicability  to  fires  where  a  strong  buoyant  field  exists 
seems  to  be  untested  as  yet.  In  previous  work,  we  used  the 
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four  flux  model  and  employed  the  gray  gas  approximation  on 
gaseous  as  well  as  liquid  pool  fires  of  moderate  size  using 
absorption  coefficients  as  input  parameters.  In  the  present 
work  they  are  computed  as  described  below,  using  the 
temperature  weighted  sum-gray  gases  model. 

The  radiative  flux  sum  F  in  the  radial  and  axial 
directions  is  given  in  the  four-flux  model  as  : 

y1  d/dy  [yFy  dFy/dy]  =  K(Fy-E)  (10) 

and 

d/dz  [T2  dFz/dz  ]  =  K  (Fz-E)  (11) 

where 

Ty  =  y/(l+Ky)  and  Tz  =  1/K  . 

Here,  K  is  the  flux  absorption  coefficient  (m_1  ). 

The  flux  absorption  coefficient  K  can  be  derived  from 
the  prevailing  concentration  and  temperature  fields  in  the 
flame,  and  the  optical  path  length  as  described  below.  The 
total  emissivity  in  the  weighted  sum-gray  gases  model  is 
evaluated  from  the  expression 

I 

ci  =  I  a e ■ , i ( T )  [1-  exp(-kiPS)]  (12) 

i  =  0 

where  the  a(i, i  represent  the  emissivity  weighting  factors 
for  the  i  ‘ h  gray  gas  based  on  the  gas  temperature,  T.  ki  is 
the  absorption  coefficient,  P  is  the  partial  pressure  of  the 
absorbing  gases,  and  S  is  the  path  length.  The  weight 
factor  a,i, i  is  physically  interpreted  as  the  fractional 
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amount  of  black  body  energy  in  the  spectral  region  where  the 
gray  gas  absorption  coefficient,  ki,  exists. 

A  relation  between  the  true  absorption  coefficient,  k, 
and  the  flux  absorption  coefficient,  K,  can  be  established 
by  assuming  thermal  equilibrium,  which  implies  that  the 
local  cell  absorptivity  and  the  emissivity  are  equal.  Then 
the  absorption  coefficient  k  is  determined  by 

k  =  cm/6n  j  (13) 

where  6n  is  the  cell  dimension.  The  flux  absorption 
coefficient  K  is  equal  to  k i  or  2k  depending  on  the 
distribution  of  the  radiation  intensity.  In  the  present 
work,  K=  2k  was  used  implying  an  isotropic  radiation 
intensity. 

4. Soot  formation  and  radiation 

In  the  absence  of  a  detailed  soot  formation  model  for 
turbulent,  diffusion  flames,  we  have  employed  the 
phenomenological  scheme  of  soot  formation  proposed  by  Tesner 
et  al,  which  was  later  applied  to  jet  flames  by  Magnussen. 
The  soot  formation  process  depends  on  the  formation  of  some 
active  radical  particles.  From  these  particles,  soot 
particles  grow.  The  formation  rate  of  radical  nuclei  is 
given  by: 

Ro,(  =  n0  +  (f-g)n  -  go  n  N  (part/m3/s)  (14) 

where  n  is  the  concentration  of  radical  nuclei,  N  is  the 
concentration  of  soot  particles,  and  n0  is  the  rate  of 
spontaneous  formation  of  radical  nuclei, 


n0  =  a0  CfU  exp  (-E/RT) 


( part/m  3/s ) . 


(  15) 


Above,  C(u  is  the  mean  concentration  of  fuel  (kg/m3),  and  f, 
g,  and  go  are  certain  constants.  The  rate  of  soot  formation 
is  expressed  by 

Rs , f  =  mP  (  a  -  bN)  n  (  Kg/m3/s)  >  (16) 

where  mP  is  the  mass  of  the  soot  particle  (Kg/part),  and  a 
and  b  are  constants. 

The  soot  and  radical  nuclei  are  formed  and  burnt  in  the 
flame  zone.  The  mean  rate  of  nucleus  and  soot  particle 
combustion  are  expressed  as  follows  : 

R„ , c  =  Rfu  n  /  CfU  (part/m3/s)  (17) 

and 

Rs.c  =  Rru  Cfi  /  Cfu  ( kg/m 3/s )  ,  (18) 

where  Cs  is  the  soot  concentration  (  Kg/m3)* 

The  absorption  coefficients  referred  to  earlier  are 
updated  for  the  calculation  of  the  overall  absorption 
coefficients  of  the  soot  +  gas  mixture  by  the  relation 

k=  k t  +  K a  -Kg  K s ,  (19) 

where  kg  is  the  absorption  coefficient  due  to  the  gas  only 
and  k<  is  the  absorption  coefficient  of  the  soot.  The  soot 
absorption  is  approximately  expressed  by: 

kg  =  4.38.  10  -«  C«  T  (20) 

5.  Solution  procedure 

The  mathematical  model  described  in  the  previous 
section  is  given  by  a  set  of  partial  differential  equations 
which  are  all  parabolic  except  for  the  axial  radiation.  All 
of  these  parabolic  expressions  are  solved  by  the  marching 
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technique  using  a  modified  GENMIX  code.  However,  the  axial 
flux  sum,  F 2 ,  is  elliptic  and  not  compatible  with  the  same 
process  and  therefore  needs  special  attention.  The  axial 
flux  sum  Fz  is  solved  separately  using  the  Tri-Diagonal- 
Matrix  Algorithm  (TDMA)  and  the  temperature  field  was 
updated.  The  GENMIX  code  employs  a  floating  grid  in  the 
radial  direction.  We  employ  50  grid  nodes  over  the  half 
width  of  the  flame.  A  maximum  axial  step  size  corresponding 
to  0.1%  of  the  current  width  of  the  flame  is  employed. 

V.  ACCOMPLISHMENTS  (  Nov.  15,  1986  -  February  28,  1988) 

During  this  contract  period  we  completed  the  following 
tasks : 

1.  Application  to  liquid  pool  fires 

Figure  1  summarizes  the  predicted  visible  flame  height 
as  a  function  of  the  reduced  calorific  output  coordinate  Q* 
for  various  liquid  pool  fires  (methanol,  heptane,  kerosene, 
and  gasoline).  Results  for  gas  flames  of  propane  and 
methane  are  also  included  here.  The  quantity  Q«  is  defined 
as 

Q*  =  [  Q  /  p  CP  Ts  g1/2  ]  2/5  (21) 

where  Q  is  the  calorific  output  of  the  flame,  p  is  the 
density  of  the  fuel  vapor,  Ts  is  the  fuel  surface 
temperature,  and  g  is  the  acceleration  due  to  gravity.  The 
visible  flame  height,  Zv,  is  obtained  from  the  centerline 
temperature  profile.  The  ratio  Zv/Q*  is  about  4.5  and  is 
consistent  with  the  value  for  buoyant  flames.  The  predicted 


Zv(m) 


gure  1  :  Comparison  of  Visible  Flame  Height  (  Zv)  Predicted 
and  Reported.  QM  is  reduced  coordinate. 


flame  heights  are  in  reasonable  agreement  with  experiment. 

We  need  to  point  out  here  that  we  did  not  calculate  the 
mass  burning  rate  for  these  liquid  pool  fires,  but  rather 
assumed  values  obtained  from  experimental  da  .a . 

2.  Determination  of  initial  boundary  conditions  for 

the  turbulent  quantities  and  application  to  propane 
flames . 

During  the  course  of  our  work,  it  was  realized  that 
optimization  of  the  initial  conditions  for  the  <-  c 
turbulence  quantities  is  needed  for  better  reproduction  of 
the  flame  behavior.  Seemingly,  the  choice  of  initial 
conditions  for  ic  and  e  are  critical  to  the  combustion  model. 
Unfortunately,  in  the  absence  of  measured  inlet  quantities, 
no  simple  or  straightforward  method  seems  to  be  available  to 
handle  this  problem.  As  a  practical  solution  to  this 
problem,  in  our  calculations  the  turbulent  quantities  were 
introduced  at  the  burner  exit  by  employing  the  following 
relations , 

no  =  cxk.Uo2,  (22) 

c0=oce.*:3/*/ri,  (23) 

where  ocr  and  a  e  are  variable  parameters,  and  n  is  the 
radius  of  the  burner  exit.  A  systematic  analysis  was 
carried  out  to  determine  those  values  for  o<g  and  ae  which 
gave  the  best  overall  flame  behavior.  Values  of  ag  in  the 
range  10-4  to  102  and  a*  in  the  range  10" 1  to  10  were 
tested.  The  final  values  chosen  in  this  study  are  those 
which  give  the  best  agreement  with  experiment  for  one  of  the 


1).  The  same  constants  were 


propane  flames  (  aK  =  10,  ae  = 
used  in  all  the  calculations  for  the  other  flames. 

Figure  2A  shows  the  variation  of  the  mean  temperature 
along  the  burner  axis  for  a  propane  flame  of  37.9  kW,  with 
and  without  radiation.  Accompanying  experimental  results 
are  shown  from  Gengembre  et  al.  The  temperature  field 
predicted  without  radiation  is  close  to  the  experimental 
results  except  in  the  peak  combustion  zone  where  it  is 
overestimated.  Also,  the  downstream  temperatures  are 
slightly  underestimated.  Upon  including  radiation,  the  peak 
temperature  as  well  as  the  downstream  temperatures  are  in 
better  agreement  with  experiment.  The  agreement  in  the 
temperature  throughout  means  that  the  combustion  zone 
height,  Zc,  and  the  visible  flame  height,  Z»,  defined  as 
above,  are  well  reproduced.  The  Wallburn  model  of  Tamanini , 
clearly  overestimates  the  temperature  field  below  100  cm. 

For  the  results  in  Fig.  2,  only  gaseous  radiation  is 
considered  and  radiation  by  soot  particles  is  not  explicitly 
included  in  the  model.  A  radiative  fraction01  XH,  is 
considered  ‘■o  account  for  the  total  heat  loss  from  the 
flame.  Incorporation  of  soot  radiation  is  considered  in  our 
further  work  on  methane  flames,  to  be  described  below. 

Figure  2B  shows  the  centerline  axial  velocity  of  the  same 
flame,  with  and  without  radiation.  The  agreement  between 
our  predicted  and  experimental  results  is  generally  good. 

3.  Inclusion  of  soot  radiation  and  application  to 


methane  flames. 


T  (K) 


Z  (cm) 

Figure  2  A:  Comparison  of  predicted  centerline  mean 

temperature  with  the  experiment  for  propane  flam 
of  37.9  kW.  Only  gas  radiation  included.  2B: 
Comparison  of  predicted  axial  velocity  along  the 
centerline  with  the  experiment 
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Figure  3A  shows  the  temperature  field  of  a  methane 
flame  along  the  centerline,  predicted  with  and  without 
radiation.  Here  soot  formation  was  considered  explicitly  by 
solving  the  differential  equations  for  soot  and  radical 
nuclei,  as  described  above.  The  absorption  coefficients  of 
the  soot  +  gas  mixture  was  calculated  using  Eq.(19).  The 
predicted  temperature  field  is  seen  to  be  close  to  the 
experimental  values  as  reported  by  Crauford  et  al .  The  soot 
concentration  field  is  shown  in  Figure  3B.  As  seen,  soot  is 
formed  in  the  fuel  rich  region  and  at  the  same  time  burnt. 
However,  in  the  present  case,  soot  is  xhausted  at  the  plume 
region,  indicating  a  smoking  tendency. 

IV.  SUMMARY 

Our  accomplishments during  this  past  contract  period 
can  be summarized  as  follows: 

1.  Application  of  the  field  model  to  liquid  pool  fires  of 
moderate  size  and  comparison  with  experimental  data. 

The  agreement  is  generally  good. 

2.  Optimization  of  -the  lc-e  initial  conditions  for  burners 
of  moderate  size,  and  model  validation  using  three 
propane  flames  of  varying  calorific  power  output.  The 
radiation  model  was  further  augmented  with  a 
temperature  weighted  gray  gas  model . 

3.  Development  of  a  model  for  soot  formation  and  soot 
radiation  and  application  of  the  model  to  methane 
flames  with  good  agreement  with  experimental  data. 


Figure  3 A: Comparison  of  centerline  mean  temperature  for 

methane  flame  of  28  kW.  Soot  radiation  included  in 
the  model.  3B:  Concentration  of  soot  (  Cg  )  along 
the  centerline  for  methane  flame. 
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These  accomp 1 i s h men t s  have  b een  communi c  a  t  e  d_  i  n 
the  following talks  and papers  : 

1.  "  Numerical  Modeling  of  Freely  Burning  Pool  Fires  ", 

paper  presented  at  the  International  Meeting  of  Fire 
Research  and  Test  Center,  Avila,  Spain,  October,  1986. 
Vol .  1,  p.  83  of  Symposium  Publication. 

2.  "  Radiative  Heat  Transfer  in  Pool-Like  Gas  Flames  ",  a 
talk  presented  at  the  Fall  Technical  Meeting  of  The 
Eastern  States  Section  of  the  Combustion  Institute,  San 
Juan,  Puerto  Rico,  December,  1986 

3.  "  Modeling  Turbulent,  Buoyant  Pool-Like  Gas  Flames 
Including  Radiative  Heat  Transfer  "  accepted  for 
publication  in  Fire  Safety  Science. 

4.  "  Numerical  Simulation  of  Pool-Like  Gas  Flames  of 
Methane  With  Special  Emphasis  on  Soot  Formation  and 
Radiation  ",  Paper  to  be  submitted  to  Fire  Safety 


Science . 


